Numerical Simulation of Mixing in Viscous - Fingering

Mixing of two fluids in viscously unstable displacements is far from being fully understood. It is not known how mixing efficiency depends on the viscosity contrast between the fluids, especially for adevction-dominated flows (Peclet number Pe > 103). It is well known that when a less viscous fluid displaces a more viscous fluid, the displacement front is unstable and leads to the formation of a pattern known as viscous fingering. However, current simulation technology is unable to cope with large viscosity contrasts (M > 30). We develop a high-resolution simulation approach that is stable for arbitrary viscosity ratios, and we study mixing under different canonical configurations with viscosity contrasts up to M = 400. We explain the observed evolution in degree of mixing through numerical simulation and dimensional analysis. We compute degree of mixing from decay in concentration variance and relate it to the stretching of material interface between the fluids due to fingering. Our analysis predicts the optimum range of viscosity contrast and Peclet number that maximize fluid-fluid interfacial area by balancing the number of fingers with their length before diffusive mixing across the sharp interface takes over. Interesting fingering patterns such as channeling and tip-splitting play an important role in this balancing act which makes degee of mixing a non-monotonic function of viscosity contrast and Peclet number. Thesis Supervisor: Ruben Juanes Title: Assistant Professor of Civil and Environmental Engineering Thesis Supervisor: Luis Cueto-Felgueroso Title: Postdoctoral Associate of Civil and Environmental Engineering 7